Tunneling "zero-bias" anomaly in the quasi-ballistic regime 
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For the first time, we study the tunneling density of states (DOS) of the interacting electron 
gas beyond the diffusive limit. A strong correction to the DOS persists even at electron energies 
exceeding the inverse transport relaxation time, which could not be expected from the well-known 
Altshuler-Aronov-Lee (AAL) theory. This correction originates from the interference between the 
electron waves scattered by an impurity and by the Friedel oscillation this impurity creates. Account 
for such processes also revises the AAL formula for the DOS in the diffusive limit. 
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A number of experiments JjJ on tunneling in disordered 
conductors consistently reveal a suppression of the differ- 
ential conductance at low biases, commonly referred to 
as zero-bias anomaly. This anomaly has found a coher- 
ent explanation in the well-known Altshuler-Aronov-Lee 
(AAL) theory of one-particle density of states. AAL 
have shown that the electron-electron interaction in the 
presence of disorder results in a negative correction to 
the density of states (DOS), which is singular at Fermi 
energy. In the case of a two-dimensional conductor, the 
AAL result reads: 
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where the parameter A > depends on the details of the 
inter-electron interaction, and is inversely proportional 
to the transport relaxation time r. Correction (Q) di- 
verges if the electron energy approaches the Fermi level, 
e — ► 0. The AAL theory assumes the diffusive motion 
of electrons, which constraints the electron energy to the 
interval e < %/t. Clearly, in the case of a strong disor- 
der, epT ~ %, this condition is not restrictive. However 
in cleaner samples (e.g., heterostructures with tunable 
density of two-dimensional electron gas Q]) the energy 
domain e > %/t becomes accessible, while the region of 
applicability of Eq. ([!]) shrinks. 

One may infer from Eq. (Q) that the correction to the 
DOS vanishes at energies larger than %/t, i.e., in the 
quasi-ballistic regime. The main purpose of this Letter 
is to show that, to the contrary, interaction does lead 
to a significant correction to the DOS even for electron 
energies exceeding the inverse transport relaxation time. 
This correction arises from the interference of scattering 
on an impurity and on the Friedel oscillation it creates. 
Account for such processes also revises the original AAL 
formula for the DOS in the diffusive limit, adding a non- 
singular however large contribution to Eq. (|l|). 

Electron density of states for energies larger than the 
inverse transport relaxation time is associated with the 
electron dynamics on the time scale shorter than r. Dur- 
ing such a small time an electron does not experience a 
large number of scattering events, i.e., the scattering on 



disorder potential can be treated in the lowest order of 
the perturbation theory, in contrast to the diffusive limit. 
This approximation accounts only for the trajectories of 
electrons which were scattered only once. We will show 
that the logarithmically divergent correction to the den- 
sity of states appears already in this approximation. 

We start with the most instructive case of the finite- 
range interaction potential, and will calculate the correc- 
tion to the one-particle DOS in the quasi-ballistic limit 
due to a single short-range scatterer. Consider an impu- 
rity at the origin; its potential U llap {v) induces a modula- 
tion of electron density around the impurity. In the Born 
approximation, one can find the oscillating correction to 
the electron density, Sn(r), which is known as the Friedel 
oscillation: 
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Here r is the distance from the impurity, v — m/Tr% 2 is 
the free-electron density of states, m is the electron mass, 
kp is the Fermi wave vector, and g = J Ui mp (r)dr. In the 
presence of interaction V(r — r') between electrons, this 
density oscillation produces an additional scattering po- 
tential, which can be presented as a sum of Hartree and 
exchange (Fock) terms: 
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where Spfri,^) ~ dn[(ri + r2)/2] is the perturbation of 
the density matrix by the impurity. The exchange inter- 
action occurs only with the electrons of the same spin, 
which is reflected by the factor 1/2 in Eq. (3b). The 
Hartree- Fock energy (|3a|) , (^) oscillates as a function of 
distance from the impurity in the same manner as Sn(r) 
does. 

The local density of states is related to the retarded 
Green function, v(e, r) = — (2/tt) ImGf (r, r). Let us find 
now the correction 5Gf (r,r) due to a coherent process, 
which includes a scattering on the impurity potential it- 
self, and a scattering on the potential (j|) formed by the 
Friedel oscillation: 
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The Green function Gf (r) for a free electron at large 
distances, k F r ^> 1, and small energies, e <C ep , is 
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in two dimensions; e is measured from the chemical po- 
tential. 

Below we will be interested in the density of states 
averaged over the spatial scales much larger than the 
Fermi wave length Xp = 2it/k F . Therefore, we should 
retain only those corrections <5^(e,r), which are smooth 
functions of r. Let us show now, using the Hartree 
contribution as an example, that Eq. (|]) indeed yields 
such a correction. This contribution corresponds to the 
following process. Electron starts motion in the point 
r, then experiences two scatterings, (1) on the impu- 
rity potential in the origin, and (2) on the potential 
formed by the Friedel oscillation in point ri, and fi- 
nally returns to point r, see Fig. |l|. Motion along this 
closed contour is represented in Eq. (Q) by the product 
Gf (r)Gf (ri)Gf (ri - r) oc exp[i^(r, n)], where 
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is the geometric phase acquired by the electron. There 
is another strongly oscillating factor in the integrand of 
Eq. (||) - the scattering potential Vh{tc\) oc sin(2k F rx). 
Obviously, the result of integration is determined by the 
domain in space where the total phase of the integrand, 
<fi(r, r\) — 2kpn, is a slow function of ri. The correspond- 
ing electron trajectories are those close to the straight 
line, see trajectory A on Fig. 1. At r\ > r, Eq. (||) yields 
the total phase of the integrand 2{e/%v F )ri. Remark- 
ably, this phase does not depend on r. As a result, the 
Hartree correction to the Green function, Eq. (|j), be- 
comes a non-oscillating function of r. Similar arguments 
can be applied to the evaluation of the exchange correc- 
tion to the Green function. The resulting expression for 
the interaction correction to the local DOS is 
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if the distance from the impurity lies within the inter- 
val m&x{\ F ,d} < r < hvp/e, drops rapidly (oc 1/r 3 ) at 
r > fivp/e, and saturates at r max{AF,<i}- Here d is 
the characteristic spatial scale of the interaction poten- 
tial, and V^(O) and V(2k F ) are the Fourier components 
of the interaction potential appearing from the exchange 
and the Hartree terms respectively. 

In order to find the averaged density of states one 
should sum up contributions of the type given by Eq. (|7]) 



from all the impurities and then average over point r 
where the correction is measured. Introducing the con- 
centration of impurities n; and using Ti/t — 2irvnig 2 , we 
arrive to the following expression for the interaction cor- 
rection to the averaged DOS in quasi-ballistic (e ^ Ti/t) 
limit: 
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with A = minjeF, Tiv F /d}. 




FIG. 1. Two typical trajectories (A, B) of an electron scat- 
tered by an impurity (/) and by the corresponding Friedel os- 
cillation (concentric arcs). The correction 5v(r) is dominated 
by the trajectories of the type A, for which the electron is 
almost scattered back at I and ri. 

In principle, the correction can have any sign depend- 
ing on the relation between V(0) and V(2kp). However, 
in any realistic system the inter-electron interaction is 
sufficiently smooth, d > X F , and 1^(0) » V(2k F ). There- 
fore, in the following we will concentrate on the exchange 
contribution, which always dominates in the correction to 
the DOS. 

The derivation of Eq. (^) is valid for energies e exceed- 
ing h/r, which is the high-energy cutoff in the AAL the- 
ory, see Eq. (0) . Quasi-ballistic formula (|^) at the bound- 
ary of the region of its applicability, e ~ Ti/t, does not 
match the AAL result ([!]). The reason of this mismatch is 
the choice of the high-energy cutoff. Physically, AAL cut- 
off means taking into account only electron states within 
a narrow energy strip, —h/r < e < 0, below the Fermi 
level, in the formation of the Friedel oscillation. This 
cut-off was dictated by the range of applicability of the 
diffusion approximation for the electron dynamics AAL 
On the other hand, our analysis leading to 
demonstrates that electron states within a much 
wider strip, —A < e < 0, are important for the correc- 
tion. As it turns out, this wider strip is important for 
the calculation of the DOS at e < Ti/t as well. To show 
this and to remedy the mismatch, below we calculate the 
DOS without using the diffusion approximation. 

We are interested in the spatially averaged density of 
states, which makes it possible to use the standard di- 
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agrammatic techniques. The correction to the averaged 
one-particle density of states has the form: 
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where e„ = n(2n + 1)T is Matsubara frequency, T is the 
temperature. (For brevity we omit the Planck constant 
in all the intermediate formulas.) We will calculate the 
correction to the electron propagator, SG(ie n ,p), in the 
first order in the screened electron-electron interaction 
V sc (i£li, Q). In the metallic regime (e ft ^> 1) the ex- 
change contribution to the propagator is: 

6G(ie n ,p) = -[G(ie n ,p)] 2 Tj2 f t|% ^(fli - <*)) 
x [T(ifl h Q)] 2 V sc (m l: Q)G(ie n - fi,, p - Q). (10) 

Here G(ie n ,p) = [ie n — £ P + (*/2 r ) signe„] _1 is the elec- 
tron Green's function in the dirty conductor, and T is 
the vertex function calculated in the ladder approxima- 
tion. As long as we are developing theory applicable for 
any relation between electron energy and H/t, we cannot 
use the usual diffusion form for vertex function. The for- 
mula valid for an arbitrary momentum Q and energy Qi 
transfer is: 
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Note, that in the limit Qi,VpQ <C r _1 , Eq. ( p"lj ) reduces 
to the standard diffusion expression. 
Integration over p in Eq. (||) gives 
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The case of a finite-range electron-electron interaction 
is especially simple, because we can replace V sc (iCli,Q) 
in Eq. ( |l2| ) by the Fourier component of the unscreened 
interaction potential V(Q). In this case the correction 
derived from Eq.(|l2|) coincides with the exchange term 
in formula (^). It means that the exchange correction to 
the DOS 
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is universal, i.e., is valid for the energies both larger and 
smaller than H/t. 

For the long-range Coulomb interaction, V(Q) — 
2ne 2 /Q, the screening should be taken into account, 
V sc (iQ h Q) = V(Q)/[l + V(Q)U(in h Q)}. Here the po- 
larization operator 
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is derived in the random phase approximation. 

Straightforward evaluation of Eq. ([l2]) with the ac- 
count for Eqs. (^lj) and ( |l4| ) yields 
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where as — I /me 2 is the Bohr radius, e = max{|e|,T}, 
and the cut-off energy A is given now by A = %Vp /as- 

Equation (|l^) gives the correction to the one-particle 
DOS due to the Coulomb interaction and is the main 
quantitative result of this Letter. The integral in Eq. (|l5| ) 
cannot be expressed in terms of elementary functions. 
However, the diffusive (e <C ^/t) and quasiballistic 
(e 3> H/t) asymptotic behaviors are easily found. 

In the diffusive limit exchange correction to the one- 
particle DOS has the form: 
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where D = v 2 f t/2 is the diffusion coefficient. The first 
term of the sum in (|l^) is the famous result of the 
Altshuler-Aronov-Lee theory [0. The second, new term 
is not singular. This part of the correction represents 
the contribution of electrons deep in the Fermi sea, with 
energies below the "%/t strip". 
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FIG. 2. The crossover function f(x), see Eq. (p.8|). 

In the quasi-ballistic limit exchange correction to the 
one-particle DOS is: 
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The leading term in the energy dependence of the cor- 
rection is oc (me) 2 at any relation between e and r. In 
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the crossover region, the correction is described by the 
formula 
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The dimcnsionlcss function f(x) is plotted in Fig. 2. 
Its asimptotes, f(x — > oo) = and f(x — * 0) = 
— (l/2)[ln(x)] 2 , allow one to obtain from ([18]) the limits 
( |j~7| ) and (|l^) respectively. 

In the practically important case of a gated semicon- 
ductor structure, the Coulomb potential is suppressed at 
all distances exceeding the separation d between the two- 
dimensional electron gas and the gate. If d 3> as, then 
the correction to the density of states coincides with Eq. 
( |l7|) at energies e > %vp I \fa~Bd, and becomes logarith- 
mical H at smaller energies, 
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with A = hvp/as- This formula is applicable for an 
arbitrary relation between e and H/t. 

The density of states Eq. (||) describes adequately the 
electron tunneling without the lateral momentum con- 
servation, such as tunneling through an inhomogeneous 
barrier. However, the electron-electron interaction affects 
the tunneling through a homogeneous barrier as well. We 
will consider below tunneling between two identical quan- 
tum wells (QW), assuming (in accordance with the ex- 
periments H) the lateral momentum conservation in the 
course of tunneling. The correction to the conductance 
we find, has a logarithmic zero-bias anomaly. 

In the absence of disorder and of the electron-electron 
interaction, the conservation of in-plane momentum im- 
plies that an electron can tunnel only if the levels of 
spatial quantization in the wells line up precisely. This 
makes differential tunneling conductance of a symmetri- 
cal double-QW system have a singular peak at zero bias. 
Disorder alone smears the singularity, leading to the I-V 
characteristic 101 
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with the width given by the inverse quantum lifetime of 
electrons in the wells, l/r s . Here V is the bias applied to 
the contact, and Go is the zero bias conductance. 

Electron-electron interaction adds a singular at zero 
bias, negative logarithmic correction to the current. De- 
tails of the calculations will be published elsewhere; the 
result for the interaction correction to the tunneling cur- 
rent is: 
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Here eV = maxjeV, T} is assumed to satisfy j8| the con- 
ditions eV <C h/r s , vp/y/aBd, and d is the separation 



between wells. Qualitatively Eq. (^l|) and Eq. (|l^) have 
the same feature: the upper cut-off for the correction is 
A = hvp/as, and by no means H/t. The oc lnel^ bias de- 
pendence, instead of esc (lneR) 2 one, appears in Eq. ( pl| ) 
because of a partial cancellation of the corrections 
coming from the intra- and inter-well electron-electron 
interaction. In the absence of interaction, Eq. ( p0| ) would 
lead to a peak in the differential conductance dl/dV at 
zero bias. Negative diverging correction (21) splits this 
peak in two. The separation between the maxima of these 
two sub-peaks is: 
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The sub-peaks should be resolved at sufficiently low tem- 
peratures, T < eV s . Estimate for V sp for the data of 
Turner et al, Ref. [0, gives V sp « 0.05mV. It is impor- 
tant, that Eqs. ( |21| ) and ( p2| ) are valid at any relation 
between eV,eV sp , and the energy H/t. 

In summary, we studied the tunneling density of states 
of interacting two-dimensional electron gas beyond the 
diffusive limit. A significant interaction-induced sup- 
pression of the density of states persists at electron en- 
ergies even larger than the inverse transport relaxation 
time, which could not be expected from the well-known 
Altshuler-Aronov-Lee theory (|]. The AAL formula for 
the density of states at low energies is also revised, and 
an additional non-singular, however large, contribution 
was found. 
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